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The Higgs field in the standard model (SM) may couple to new physics sectors related with dark 
matter and/or massive neutrinos. In this paper we propose a novel signature, the boosted di-Higgs 
boson plus (Kt (which is either a dark matter or neutrino), to probe those new physics sectors. In 
a large class of models, in particular the supersymmetric SMs and low scale seesaw mechanisms, 
this signature can play a key role. The signature has clear background, and at the ^/s =14 TeV 
high luminosity (HL-)LHC, we can probe it with production rate as low as ~ 0.1 fb. We apply it 
to benchmark models, supersymmetry in the bino-Higgsino limit, the canonical seesaw model and 
little Higgs model, finding that masses of Higgsino, right-handed neutrino and heavy vector boson 
can be probed up to ~ 500 GeV, 650 GeV and 900 GeV, respectively. 

PACS numbers: 


New physics below the iceberg Discovery of Higgs 
boson at the large hadronic collider (LHC) [I[ completes 
the standard model (SM). At the same time, however, it 
opens a new era for particle physics: this new resonance 
might be just a small tip of a big iceberg, and below it 
could hide a mystery of new world. Looking around this 
small tip, we may find clues for new physics. Actually, 
we do have convincing arguments to support this belief 
from several motivations for new physics beyond the SM 
(BSM). 

The most common theoretical argument for new 
physics is from the notorious gauge hierarchy problem 
caused by the quadratic divergence of (Higgs mass)^ pa¬ 
rameter. Solutions to this problem introduce new par¬ 
ticles coupled to the SM Higgs field. For instance, in 
the supersymmetric SMs (SSMs) the Higgs doublet 
(along with Hd by virtue of anomaly cancelation) partic¬ 
ipate in quite a few interactions. The second argument 
for new physics is from dark matter (DM), whose interac¬ 
tions with the visible sector may be through the SM Higgs 
boson h (Higgs portal DM). Last but not the least, neu¬ 
trinos in the renormalizable SM are massless in conflict 
with the observation. Mechanisms, such as various see¬ 
saw mechanisms, to generate neutrino masses may again 
introduce new couplings to h. It is well expected that 
these new interactions of h could contribute to the pro¬ 
duction of h (plus something else) at LHC. In this letter 
we will concentrate on the di-Higgs production, which is 
common in these mentioned BSM contexts, has a bright 
prospect at LHC and moreover provides a new angle on 
di-Higgs physics. 

BSM with boosted di-Higgs bosons plus Di- 

Higgs boson search is one of the focuses in the upcoming 
LHC run, aiming at examining the Higgs potential 
On the other hand, on top of the di-Higgs, a large class 
of BSM models can produce associated objects X, e.g., 
A is a missing particle at LHC such as DM or neutrino. 


Moreover, both Higgs bosons may be boosted if they are 
produced from heavy particle decay. These additional 
features could greatly enhance the di-Higgs searches at 
a hadronic collider. In the following we select two well- 
known BSM models where A is a DM and neutrino, re¬ 
spectively. 

The first example is the neutralino sector of the mini¬ 
mal SSM (MSSM) with bino-Higgsino accessible at LHC 
only, ijj = {B, H°, H^). ^ Since we are interested in 
boosted di-Higgs, we assume such a mass hierarchy: Hig- 
gsinos ^ (with mass /i ^ 500 GeV) are much heavier 
than bino B (with mass < 100 GeV). The mass eigen¬ 
states Xi, 2 , 3 , respectively having masses Mi_ 2,3 in ascend¬ 
ing order, are related to ip via ipi = ZpjPxj- X 2 and xs 
are Higgsino-like, constituting a pseudo Dirac fermion 
pair with masse splitting of a few GeVs; xi is treated as 
massless for the moment. Electroweak gauge interactions 
of these particles are described by 

^ruXza'PjXjh + ( 1 ) 

^ ^ C/CJo ^uj 

with 9yj the Weinberg angle and oT = 

Z^Z^ — Z'^Z'^. We have assumed CP-invariance and an 
exact decoupling between two Higgs doublets. Note that 
the diagonal couplings al^ are suppressed by the small 
mass splitting between M 2 and M 3 . 

The interesting signature is produced along the chain 
pp ^ Z* ^ X 2 (-t Xi^)X 3 (-t Xi^)- In the high energy 
limit MI 3 — Mf TO^, we can apply the Goldstone 
equivalence theorem and obtain the following equations 
for decay branching ratios of X 2,3 (We refer to [3] for 


^ The gravitino-Higgsino system with very light (says '^keV) grav- 
itino from the low scale supersymmetry-breaking models fit this 
scenario better. We leave it for a future study. 
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relevant discussions on this point): 

r(x2 fe + xi) _ r(x2 -» z + xi) 
r(x3 ^ + xi)"" r(x3 -^h + xi)' 

This theorem will be used frequently in this paper. In 
order to approach the maximal branching ratio of the di- 
Higgs boson channel, r(xz h+xi) ~ r(xi ^+xi) = 
50% is favored. And it holds when the Higgs sector gives 
a large tan/3 limit and moreover |^| k, |M 2 , 3 | /§> Mi. In 
this limit the cross section of di-Higgs plus can reach 
1.2 fb for 500 GeV Higgsinos. 

The second example is the heavy right-handed neu¬ 
trino (RHN) N in the type-I seesaw mechanism imple¬ 
mented in the local B — L model Q . The minimal model 
is described by the following Lagrangian: 

—^N =gB-LNXfj,PiiNZ^_^ + yMiHPfiN 

+ ^XN‘i>N^N + h.c. + V{^,H), (3) 

where we consider only one family of RHN for simplicity. 
N carries one unit oi B — L charge and gains Majorana 
mass from the coupling to the B — L Higgs field 3>, which 
develops VEV from the potential V ($, H) and breaks 
U{1)b-l attheTeV scale. The term inR($,i?) 

mixes $ and the neutral component of SM doublet 
H, leading to two Higgs bosons in the mass eigensate: 

h = cos 9 H'^ +sin 6 ^ 11 , (j) = — sinO + cosO^b.. (4) 

We have decomposed $ = {v^f, + + i4>/)/-\/2 similarly 

for H^. The current data tells that h is quite SM-like, 
and the measured Higgs signal strength imposes an up^r 
bound on the mixng angle sin0 < 0.34 at the 95% CL [Gj. 
For such a small mixing angle, the SM-like Higgs boson 
mass squared can be approximated as 

ml Ri Xhv'^ - sin^ 6»m^, (5) 


where Xh is the usual Higgs quartic self coupling. Let 
us mention that the heavier <f) and the larger 9 are good 
for solving the metastability problem of Higgs potential 
near the Plank scale, which is attributed to the relatively 
small Xh = Asm ~ 0.26. For example, for = I.O TeV 
and sin0 = 0.3, we now need Xh ^ 1-7. Note that such 
discussions can be easily generalized to other models. 

The heavy N can be abundantly pair-produced via 
the resonance Zb-l/4>'- PP Zb-lI4’ NN. The 
Zb-l channel is suppressed by the small branching ra¬ 
tio Bt{Zb-l —>■ NN) K. 3%. In addition to that, it is 
strongly constrained by the di-lepton resonance search 
PP Zb-l at LHC Q - The (j) channel is due to 
the gluon-gluon fusion production of </), described by the 
usual dimension-five operator for (j) —>■ gg: 


jiv T 

A 994’'-'uu'^ I 


( 6 ) 


with Cgg(f, = — sin0Q;s/27r. However, this effective de¬ 
scription becomes invalid when is much heavier than 
the top quark mass. Therefore, in the actual LHC anal¬ 
ysis, we will utilize the results from the CERN Yel¬ 
low Report Q which takes into account the finite top 
quark mass effect. For Aat ^ 1 one can naturally expect 
Br(^ NN) ~ 100%. The RHN pair is followed by the 
decay N ^ vlH with vl the active neutrino, giving rise 
to the boosted di-Higgs plus signature. Because RHN 
is heavy, the branching ratios satisfy the relations ii: 

T{N ^ i^Lh) ^ r(7V ^ ulZ) = ir(/V ^ iW), (7) 

which again follows from the equivalence theorem. For 
Mat = 0.5 TeV, one can reach a signal cross section 1.5 
fb given that m^j, is about 1 TeV. 

Comments are in order. The conventional single RHN 
(below the weak scale) production counts on sizable mix¬ 
ings between the sterile and active neutrinos [3 El, 
which would be extremely small in generic cases. But 
probably RHN has gauge and/or Yukawa couplings, e.g., 
in the local B — L models, and thus the pair production 
of heavy RHN is hopeful at 14 TeV LHC or other future 
colliders [l^. In particular, to our knowledge, using a 
scalar resonance is novel in RHN production and may of¬ 
fer the unique chance to probe type-I seesaw mechanism. 
Simplified models for 2h+fiT There are many other 
well motivated BSM models that can produce this sig¬ 
nature; see an incomplete collection in the supplement 
material. Therefore, before we head toward the detailed 
collider study, this would be the right place to develop 
simplified models which could be used for a more general 
study of this signature. 

Two classes of models are of interest. Let us start 
with the first class, where the missing particle is DM, or 
more widely a neutral particle stable at the collider time 
scale. The dark sector is supposed to consist of several 
dark states; we will use this term to genetically refer to 
particles undiscovered. And we will consider two dark 
states for concreteness. Then their interactions with the 
Higgs boson would take the following forms: 

A//iXiX2, PshSiS2, gzhZ^{Zi)g. (8) 

(Xi, Si, Zi) is a (Majorana fermion, real scalar, vector 
boson) DM with negligible mass; the heavy dark state 
X 2 , etc., is in the sub-TeV mass region. Similarly, the 
second class with neutrino as the missing particle can be 
built. RHN-like heavy state N couple to neutrino and 
Higgs boson via the effective operator 


XjyhNvB- (9) 

With these, the event topology boosted 2h-\-^T is gener¬ 
ated through pair production of heavy dark states which 
decay into DM/neutrino plus Higgs bosons. 


V 
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In the simplified models we do not specify the produc¬ 
tion mechanism for the father particle F such as X 2 since 
it is fairly model dependent. We can consider three pos¬ 
sibilities for F productions. Firstly, the father particle F 
itself such as the Higgsino has electroweak interactions so 
that the F pair can be produced via Drell-Yan processes. 
This mechanism does not involve extra particles and cou¬ 
plings. Secondly, in the presence of a new particle Y cou¬ 
pling to light quarks q as ^ qYF (just schematically), F 
can be pair-produced via exchanging Y in the t—channel. 
In the supplement, the vector F pair production in the 
little Higgs model is discussed in this way. Last but not 
the least, if Y has interactions like ^ Yqq/GG -I- YFF, 
then the F pair can be resonantly produced. But as 
mentioned before, this contribution will be stringently 
restricted if Y also couples significantly to leptons. The 
latter two possibilities may provide more effective pro¬ 
duction mechanisms than the first one, provided that the 
mass and couplings of Y (to quarks and F) are proper. 
In the following studies, a proper production mechanism 
for the F pair will be assumed. 


Boosted 2h -|- at the 14 TeV LHC The hadronic 
46 mode is dominant in the di-Higgs decay, and becomes 
relieved from the huge QCD backgrounds (BGs) in the 
boosted region. For example, D. Lima et al. show that 
with the help of the boosted di-Higgs channel one can 
finally reach 1.2 x Asm at 14 TeV HL-LHC 1^ . ^ This 
channel was believed to be not promising 13, Hj outside 
the boosted Higgs region. In some new physics scenar¬ 
ios, where the di-Higgs is produced from a heavy reso¬ 
nance decay, this channel is even more remarkable [ 13 - 
[l9l| . Other related studies involving boosted (maybe ex¬ 
tra) Higgs boson(s) are also motivated in various con¬ 
texts [isl. [2iil - E3 |. In our signature, the (highly) boosted 
di-Higgs is further strengthened by a large missing en¬ 
ergy, that could lead to the much earlier discovery of 
2 h+^T than the previous signatures without fir ■ ^ 


The BGs of 2h+fiT are similar to those of the di-Higgs 
signature: the irreducible BGs (dominated by QGD 46 
and Zbb) and the reducible BGs (dominated by the semi- 
leptonic tt pair); 4j with j = g,u,d,s,c can be mis- 
tagged as 46, but the cross section is < 10“^ x (t( 46) 
and thus negligible. In the irreducible BGs, is due to 
the limited detector resolution, and the distribution of 
depends on the detector setup: e.g., at 7 TeV LHC 
it respects a Gaussian distribution with the central value 
~ 0.5i/^ Pt GeV [i^. While in the semi-leptonic tt BG, 
light flavor jets can be mis-tagged as 6-jets and at the 


^ The ATLAS search on data set of 8 TeV 19.5 ffa—^ shows that 
the current search sensitivity to the boosted di-Higgs channel is 
around an order of magnitude above the SM di-Higgs production 
rate [T^ . without using the substructure technique. 

® CMS searched for this signature without using boosted Higgs- 
tagging fSSl . They required that all 4b’s have to be resolved. 


same time the leptons may be missed, which happens, 
besides owing to the limited lepton tagging efficiencies 
(especially for r), if the leptons go outside the kinematic 
region, i.e., pt{1) > 10 GeV and \r]\ < 2.5, or they are not 
isolated, i.e., the scalar px sum of particles in the vicin¬ 
ity (Ai? < 0.5) of the lepton is greater than 10% x pt{1)- 
MadGraph5_aMC@NLO [1^ is used to calculate the LO 
cross sections for 46 (861 pb) and Zbb (1 09 p b at NLO 
after multiplying a constant iC-factor 1.4 [23]); the bot¬ 
tom quarks in these BGs are required to have pt(6) > 20 
GeV, |? 7 ( 6 )| < 2.5 and Ai?( 6 , 6 ) > 0.4. The NNLO cross 
section of the semi-leptonic tt including all lepton flavors 
can be found in Ref. [^, 382 pb; no cuts are imposed on 
the top decay products. 

We take heavy RHN pair production via the 
(/)—resonance as the representative signal process, set¬ 
ting m^p = 2Mjv -b 30 GeV and Br(A —>• hvx) = 1. 
FeynRules [1^ is used for generating the UFO model 
files of the model [s^. The signal (and BGs) events 
are generated by MadGraph5_aMC@NLO, and passed to 
PythiaG [33 for SM particles decay, parton showering 
and hadronization. The detector effects are simulated 
with Delphes3 [s^j, where we choose the default ATLAS 
detector setup. The 6 -tagging efficiency is set to 0.7 with 
mis- 6 -tagging rates for c- and light-flavored jets assumed 
to be 0.1 and 0.015 [s^, respectively. Then, the par¬ 
ticle flow information from Delphes3 is analyzed using 
fastjet 341. 

Gomments on the generalization of the above represen¬ 
tative model are in orders. Firstly, the invisible particle is 
assumed to be massless. However, our results can be gen¬ 
eralized to the massive invisible particle case {fix refers 
to the measurable ptx in this case), because the variables 
in our analysis actually are merely sensitive to the mass 
difference between the invisible particle and the father 
particle. Secondly, we have chosen the Yukawa coupling 
of (j) such that the narrow width approximation holds 
(F^ < O.Im^); moreover, mp — 2 Mn is small. So, it is 
justified to conclude that the kinematic properties of the 
RHN pair produced here are similar to those of the F 
pair which is produced in other ways through an off-shell 
s-channel mediator or a t-channel mediator. 

In selecting the signal events, lepton veto is imposed 
first. Next, the Gambridge/Aachen algorithm 33 with 
an appropriate cone size i?fat is used to recluster fat- 
jets. Then the BDRS algorithm [s^ is applied to resolve 
their substructure. More concretely, a Higgs jet can¬ 
didate during the declustering should have large mass 
drop fi = < 0.67, and not too asymmetric mass 

splitting y = > q.OO. After filter¬ 

ing, the three hardest sub jets inside the fat-jet, which 
are reconstructed by the anti-kt algorithm with Rfm = 
min(0.3, i?{, 5 / 2 ), are identified as the ingredients of the 
Higgs jet candidate; also the two leading subjets are re¬ 
quired to be 6 -tagged. 
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The separation angle between b and b from Higgs de¬ 
cay, namely the Higgs jet cone size, depends on the en¬ 
ergy of Higgs boson. For each RHN mass Mjv, we scan 
^fat € [0.5,3] with step size 0.1 and find that 

retains most of the signal events whose leading two fat- 
jets pass the BDRS Higgs-tagging criteria and the cut 
™Higgs S [85,160] GeV (A more refined mass interval for 
a given Mjv is given in Table IL). In Table |T] we show the 
options of for different values of Mn- The fraction 

of retained signal events (selection efficiency) and cross 
sections of BGs after Higgs-tagging are also displayed. 
As expected, a smaller 1?^““ should be adopted for a 
heavier N. It then greatly reduces the probability that 
two close QGD jets mimic the Higgs jet, thus suppressing 
BGs (especially the QGD BG from 46’s). 


Mjv (GeV) 

omax 

-^fat 

S 

tt (fb) 

%66(fb) 

QCD(fb) 

200 

2.0 

2.1% 

5.2 

13.9 

591.8 

300 

1.8 

2.5% 

3.9 

9.2 

399.4 

400 

1.6 

3.2% 

2.4 

5.6 

241.7 

500 

1.6 

3.9% 

2.4 

5.6 

241.7 

700 

1.4 

5.4% 

1.3 

2.9 

117.1 

1000 

1.0 

6.8% 

0.46 

0.41 

14.5 

2000 

0.8 

8.8% 

0.19 

0.1 

3.2 


TABLE I: The selection efficiency S and cross sections of BGs 
after applying Higgs-tagging which takes mniggs € [85, 160] 
GeV and = A^r- 

There are some other powerful discriminators in our 
analysis, such as fir and the reconstructed two Higgs 
boson masses. Requiring fJx ^ 100 GeV after Higgs- 
tagging leaves the semi-leptonic tt as the dominant BG. 
In addition, since we consider pair production of heavy 
particles, the Mt^ variable [s^, which reflects mass differ¬ 
ence between the father particle and its invisible daugh¬ 
ter, is quite useful. According to these variables, we de¬ 
sign four signal regions (SR) which are suitable for differ¬ 
ent Mjv; see Table m where the name of SR is indicated 
by the RHN mass. The cuts of each SR are optimized 
with respect to the corresponding benchmark point. Note 
that SR500 will be used for Mjv above 500 GeV because 
its BGs already become negligible. 


signal region 

SR200 

SR300 

SR400 

SR500 

Selection cuts 

Lepton veto & tan veto 

Two Higgs tagged jets 

$j./GeV 

> 100 

> 100 

> 200 

> 300 

mnJGeV 

[90,150] 

[90,150] 

[90,150] 

[90,150] 

mjij/GeV 

[80,140] 

[70,150] 

[80,140] 

[80,140] 

Mt2 {hi, h2) 

> 130 

> 170 

> 215 

> 300 


TABLE II: Optimized cuts for the four signal regions. 


How far can we reach? For a given RHN mass, the 
signal cross section that provides at least 3-cr significance, 
i.e., S/y/T+B > 3, is 

1 9 -l-1/81 -l- 36 X . 

^ iir-5-' 

where C is the integrated luminosity and a-re cuts 
efficiencies for the signal and BGs, respectively. In Fig.[T] 
we show the reach limit of our representative signal pro¬ 
cess in the light green shaded region. One can see that 
the production rate as low as 0(0.1) fb can be reached 
for Mjv > 500 GeV. 


C = 3000 fb ^ Reach 



M/(GeV) 

FIG. 1: Based on the representative signal process, the shaded 
region can be probed with 3-a significance at 14 TeV 3000 
fb~^ LHC. We display the cross sections of di-Higgs-|-.^T in 
four models: (I) seesaw (orange), (II) MSSM (blue), (III) su¬ 
persymmetric seesaw (cyan) and (IV) little Higgs (red) mod¬ 
els. 

Now we explain the above model-independent reach in 
several benchmark models, where the kinematic proper¬ 
ties of the father particle (F) pair are similar to those 
of the RHN pair in the representative signal process. 
They are: (I) the neutrino system in the Type-I see¬ 
saw that presents a scalar resonance </> with a large mass 
> 2Mjv (For such heavy (j), its production rate is 
obtained from Ref. Q rather than Eq. ®.), the mixing 
angle sin0 = 0.3 and Br(^ —>■ TV TV) = I; (H) the bino- 
Higgsino system in MSSM with Br(iV° ^ Bh) = 50%; 

(HI) the sneutrino system vg — N'm the supersymmetric 
seesaw with the left-handed sneutrino vg produced 
via the Drell-Yan process and Br(?x N + h) = 50%; 
(IV) the massive vector boson system Zh — Ah in the 
little Higgs model m, IS] with Zh produced by the 
t-channel exchange of a quark partner and Br(Z// —>■ 


^ Converting the signal reach to the limit on the mixing angle for 
we get sin0 > 0.23 for around 800 GeV. But it is can be 
improved much at the 100 TeV collider [^ . 



























5 


Ah + h) = 100%. We put the details of the later two 
models in the supplement. The S-tr discovery lines are la¬ 
belled in Fig. [H We observe that (RHN, Higgsinos, Zh) 
with mass up to ~ (0.50, 0.65, 0.9) TeV can be probed 
at 3-tT level, whereas the sneutrino case can be hardly 
probed. Note that the invisible particle is assumed to be 
massless, and the violation of this assumption will lead to 
a shift of these curves toward left, roughly by the invisible 
particle mass. 

Conclusion In the light of a large class of BSM mod¬ 
els, we propose a noble signature, namely the boosted 
di-Higgs plus to probe new physics. Good prospects 
of models like MSSM, type-I seesaw and little Higgs are 
demonstrated, finding that masses of Higgsino, right- 
handed neutrino and heavy vector boson can be probed 
up to ^ 500 GeV, 650 GeV and 900 GeV, respectively. 
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